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ABSTRACT 

Dynamical Dark Energy (DE) is a viable alternative to the cosmological constant. 
Yet, constructing tests to discriminate between A and dynamical DE models is difficult 
because the differences are not large. In this paper we explore tests based on the galaxy 
mass function, the void probability function (VPF), and the number of galaxy clusters. 
At high z the number density of clusters shows large differences between DE models, 
but geometrical factors reduce the differences substantially. We find that detecting a 
model dependence in the cluster redshift distribution is a hard challenge. We show that 
the galaxy redshift distribution is potentially a more sensitive characteristics. We do 
so by populating dark matter halos in N— body simulations with galaxies using well- 
tested Halo Occupation Distribution (HOD). We also estimate the Void Probability 
Function and find that, in samples with the same angular surface density of galaxies 
in different models, the VPF is almost model independent and cannot be used as 
a test for DE. Once again, geometry and cosmic evolution compensate each other. 
By comparing VPF's for samples with fixed galaxy mass limits, we find measurable 
differences. 
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1 INTRODUCTION 

High redshift supernovae, anisotropies of the cosmic mi- 
crowave background (CMB), as well as data on the large- 
scale galactic distribution (Riess et al. 1988, Perlmutter et al 
1988, Tegmark et al. 2001, De Bernardis et al 2000, Hanany 
et al 2000, Halverson et al 2001, Spergel et al 2003, Percival 
et al. 2002, Efstathiou et al 2002) indicate that ~ 70 % of the 
world contents are due to a smooth component with largely 
negative pressure. This component is dubbed dark energy 
(DE). Recently Maccio, Governato & Horellou (2005) pre- 
sented further arguments in favor of DE based on the lo- 
cal (~ 5Mpc) Hubble flow of galaxies. The nature of DE 
is still open for debate. Candidates range from a positive 
cosmological constant A - yielding a ACDM cosmology - to 
models with a slowly evolving self-interacting scalar field (j> 
(dynamical DE; Ratra & Peebles 1988; Wetterich 1988) to 
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even more exotic physics of extra dimensions (e.g., Dvali & 
Turner 2003). 

ACDM cosmologies are easy to study and fit most data. 
Unfortunately, to give a physical motivation to the value of 
A, we need a fine-tuning of vacuum energy at the end of the 
last phase transition. To rival the success of ACDM other DE 
models ought to predict observables hardly distinguishable 
from it, so that discriminatory tests on DE nature are not 
easy to devise. 

Up to now, most tests based on large scale structure 
dealt with the evolution of the cluster distribution. Using 
the Press-Schechter-type approximations (Press & Schechter 
1974, PS hereafter; Sheth & Tormen, 1999, 2002 ST; Jenk- 
ins 2001) the expected dependence of the cluster mass func- 
tion on DE nature was extensively studied (see, e.g., Wang 
& Steinhardt 1998, Haiman et al 2001, Majundar & Mohr 
2003, Mainini, Maccio & Bonometto 2003). The results were 
used to predict the redshift dependence of various observ- 
ables, such as temperatures (T) or photon counts (N). In 
Sec. 2, we compare the (virial) mass functions for different 



© 2005 RAS 



2 Solevi et al. 



DE. An important - and often overlooked - factor is the 
dependence of halo concentration c on the DE equation of 
state. For given virial mass the concentration varies with DE 
nature up to 80 %, as shown by simulations (Klypin, Maccio, 
Mainini & Bonometto 2003; Linder & Jenkins 2003; Kuhlen 
et al. 2005). The model dependence of c is so strong that 
it can be used as a possible discriminatory test for strong 
lensing measurements (Dolag et al 2004, Maccio 2005). 

In this paper we shall focus on galactic ~ 10 12 /i _1 Mq 
scales which are also interesting for testing models of DE. 
In order to make a prediction we first need to know how 
to generate a distribution of galaxies, not just dark matter 
halos. There are different ways for doing this. We decided 
to use recent results on the Halo Occupation Distribution 
(HOD): a probability to find N galaxies in a halo of mass 
M. The HOD properties have been studied in details (Sel- 
jak 2000; Benson 2001; Bullock, Wechsler & Somerville 2002; 
Zheng et al. 2002; Berlind & Weinberg 2002; Berlind et al. 
2003; Magliocchetti & Porciani 2003, Yang, van den Bosch 
& Mo 2003; van den Bosch, Yang & Mo 2003). Numeri- 
cal simulations including gas dynamics (White, Hernquist 
& Springel 2001; Yoshikawa et al., 2001; Pearce et al. 2001; 
Nagamine et al. 2001;Berlind et al. 2003, Yang et al. 2004) 
or semi-analytical models of galaxy formation (Kauffmann, 
Nusser & Steinmetz 1997; Governato et al. 1998; Kauffmann 
et al. 1999a,b; Benson et al. 2000a,b; Sheth & Diaferio 2001; 
Somerville et al. 2001; Wechsler et al. 2001; Benson et al. 
2003a; Berlind et al. 2003) were used, to find a law telling 
us how halos split in sub-halos hosting individual galaxies. 

In this way we formulate predictions on galaxy mass 
functions and their z-dependence. In order to compare pre- 
dictions with data we then need to disentangle the evolution 
of the mass function from the evolution of the M/L ratio, be- 
cause observations provide luminosities of galaxies and our 
estimates give us halo masses. In this respect, one of the aims 
of this work is to estimate how precisely the Mj L evolution 
should be known, in order to use data on galactic scales to 
test DE nature. Potentially, M/L evolution can be predicted 
using galactic evolution models (see, e.g., Bressan, Chiosi & 
Fagotto 1994, Portinari, Sommer-Larsen & Tantalo 2004, 
and references therein). Such predictions can be compared 
with weak lensing results or satellite dynamics (Prada et al 
2003). The latter methods will provide estimates of virial 
M/L for samples of galaxies at different z. This is exactly 
what we need in order to test different models for DE. Here 
we find that, for some statistics, the expected signal, i.e., 
the differences between models are rather large. So, there is 
a hope to detect DE effects in spite of uncertainties in M/L 
ratios. Certainly, if galactic evolution predictions and high-z 
M/L estimates are compared, one can hardly prescind from 
taking into account accurately the impact of DE nature. 

We run a series of N— body simulations with different 
equations of state. In these models the ratio w = Pde/Pde of 
the DE pressure to the energy density varies with z accord- 
ing to field dynamics. The models considered here are the 
Ratra-Peebles (RP, Ratra & Peebles 1988) and the mod- 
els with the supergravity (SUGRA, Brax & Martin 1999, 
2001; Brax, Martin & Riazuelo, 2000). Appendix A gives 
a short summary on these models. Each model is speci- 
fied by an additional parameter - the energy scale A of the 
self-interacting potential of the scalar field. Here we take 



A = 10 3 GeV for both models. In RP (SUGRA), w shows 
slow (fast) variations with z. 

The void probability function is an obvious candidate 
for discriminating models. Fluctuations grow differently in 
the models. So, one may expect some differences in VPF. 
We use galaxy distributions to estimate the void probability 
function (VPF) at different redshifts. While measuring VPF 
in simulations is straightforward, mimicking observations is 
slightly more complicated because it requires corrections for 
geometrical effects and because the answer depends on a 
definition of galactic populations. At z = no model de- 
pendence of the VPF is expected or found. Predictions at 
higher z depend on how galaxy samples are defined. In par- 
ticular, we show that, if equal angular density samples are 
considered, VPF results are independent of the DE nature. 
On the contrary, if we select samples above fixed galactic 
mass M, a significant signal is found, which can be useful 
for testing the DE nature. 

The plan of the paper is as follows: In Sec. 2 we discuss 
how geometry and galactic evolution affect the redshift dis- 
tributions of galaxies and clusters. In Sec. 3 and 4 we discuss 
the simulations and prescriptions of populating halos with 
galaxies. In Sec. 5 and 6 results on redshift distributions and 
the VPF are given. Sec. 7 is finally devoted to discussing our 
results and future perspectives. 



2 GEOMETRICAL AND EVOLUTIONARY 
FACTORS 

Let us consider a set of objects whose mass function is n(> 
M,z). In a spatially flat geometry, their number between z 
and 2 + Az, in a unit solid angle, is given by 

pz+Az 

N(> M, z, Az) = / dz' D(z') r 2 (z') n(> M, z) (1) 

J Z 

with D(z) = dr/dz. For flat models: 

Here Q m (z) is the matter density parameter at the redshift 
z. The Friedman equation can be written in the form 

rr1,s 87T^ Pmo (l + Z ) 3 ^ 2 O mo (l + zf 

H ^ = Y G n m (z) =H ° n m (z) ■ (3) 

Along the past light cone, adr = —cdt. So, by dividing 
the two sides by dz = —da /a 2 , one finds that a dr/dz = 
a 2 cdt /da. Accordingly, D(z) — dr/dz — c/H(z) is derived 
from eq. (3). 

When w is a constant, a useful expression 

D 2 (z) = (c/H )(l + z)- 3 [Q mo + (1 - fi mo )(l + zf™]- 1 (4) 

can be obtained, which allows one to see that the geometrical 
factors increase both when to decreases and fl mo decreases 
in models with w < 0. 

An extension to dynamical DE can be performed either 
by using the interpolating expressions yielding £l m (z), for 
RP and SUGRA models, provided by Mainini et al (2003b) 
or, equivalently, through direct numerical integration. We 
obtain the geometrical factor r 2 (z)D(z) shown in the upper 
panel of Fig. 1. 
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Figure 1. Geometrical and evolutionary terms on the cluster 
mass scale. In both panels ACDM with f2 mo = 0.2 is above 
ACDM with n mo = 0.3. On the contrary, in the upper and lower 
panels RP and SUGRA lay on the opposite sides of ACDM. In 
the latter case, cancellation between geometrical and evolutionary 
effects is therefore expected. 
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Figure 2. The upper panel shows how the redshift distribution 
on the cluster mass scale depends on the models. The cancel- 
lation expected from Fig. 1 has occurred and the models with 
!l mo = 0.3 are quite close, while ACDM grantes large differences 
when varying Q m o- In the lower panel, similar curves are shown 
for halos on galactic mass scales. In this case, evolution is almost 
model independent and the geometrical factor causes an apprecia- 
ble difference. On these scales a lower Q m o yields a different halo 
density, unless the spectrum is normalized to unphysical levels. 
Hence, only models with Q m o = 0.3 are shown. 



is 70km/s/Mpc in all models (h = 0.7). In the absence of 
number density evolution, the upper panel of Fig. 1 shows 
the dependence of the angular number density on z. 

In the PS formulation, the expected differential cluster 
number density n(M), at a given time, is then given by the 
expression 

M 

f{v)vd\ogv= — n{M)MdlogM . (5) 
p m 

Here p m is the matter density, v = 5 c /um is the bias factor, 
M is the mass scale considered, om is the r.m.s. density 
fluctuation on the scale M and 8 C is the amplitude that, in 
the linear theory, fluctuations have in order that, assuming 
spherical evolution, full recollapse is attained exactly at the 
time considered (in a standard CDM model this value is 
~1.68; the difference, in other model, ranges around a few 
percent). As usual, we took a Gaussian f(v) distribution. 

Together with eq. (5), we must take into account the 
virialization condition, which yields significantly different 
density contrasts A„ in different DE models. Further details 
can be found in Mainini et al (2003a). 

In the lower panel of Fig. 1 we show the evolution of 
the number of halos of mass > W 14 h~ 1 MQ, in comoving vol- 
umes. All models are normalized to the same cluster number 
today and the redshift dependence of n(> M,z) is clearly 
understandable, on qualitative bases: When ACDM models 
are considered, the evolution is faster as we approach stan- 
dard CDM. On the contrary, RP and SUGRA yield a slower 
evolution than ACDM. 

The important issue is that, while both the geometrical 
factor and the evolutionary factor of ACDM (with f2 mo = 
0.2) lay above ACDM (with Q, mo = 0.3), RP and SUGRA 
factors lay on the opposite sides of ACDM. 

When the two factors are put together this causes the 
effect shown in the upper panel of Fig. 2, a strong signal 
on Q mo and a widespread cancellation for DE models, com- 
pared to ACDM. Discriminating between different DE na- 
tures, from this starting point, is unavoidably a hard chal- 
lenge. 

Geometrical factors do not depend on the mass scale. 
Instead, evolutionary factors are known to have a stronger 
dependence on the model for larger masses. As cancellation 
is almost complete on cluster scales, it is to be expected that 
geometrical factors yield a significant signal on lower mass 
scales. This is shown on the lower panel of Fig. 2, where 
halos of 10 12 h^M® are considered. A halos of this mass is 
expected to host a normal galaxy. More massive halos are 
expected to host many galaxies. Hence, this plot cannot be 
directly compared with observations. Its main significance is 
that such lower mass scales deserve to be inspected because 
DE signals are expected to be strong enough on these scales. 



3 SIMULATIONS 

The simulations run for this work are based on a ACDM 
model and two dynamical DE models, with the same mat- 
ter density and Hubble parameters (fi mo = 0.3 and h = 0.7). 
The simulations are run using the Adaptive Refinement Tree 
code (ART; Kravtsov, Klypin & Khokhlov 1997). The ART 
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code starts with a uniform grid, which covers the whole com- 
putational box. This grid defines the lowest (zeroth) level of 
resolution of the simulation. The standard Particles-Mesh 
algorithms are used to compute density and gravitational 
potential on the zeroth-level mesh. The ART code reaches 
high force resolution by refining all high density regions us- 
ing an automated refinement algorithm. The refinements are 
recursive: the refined regions can also be refined, each subse- 
quent refinement having half of the previous level's cell size. 
This creates a hierarchy of refinement meshes of different 
resolution, size, and geometry covering regions of interest. 
Because each individual cubic cell can be refined, the shape 
of the refinement mesh can be arbitrary and match effec- 
tively the geometry of the region of interest. 

The criterion for refinement is the local density of par- 
ticles: if the number of particles in a mesh cell (as estimated 
by the Cloud-In-Cell method) exceeds the level n t hrosh, the 
cell is split ("refined") into 8 cells of the next refinement 
level. The refinement threshold may depend on the refine- 
ment level. The code uses the expansion parameter a as 
the time variable. During the integration, spatial refine- 
ment is accompanied by temporal refinement. Namely, each 
level of refinement, I, is integrated with its own time step 
Aci; = Aoo/2', where Aao is the global time step of the 
zeroth refinement level. This variable time stepping is very 
important for accuracy of the results. As the force resolution 
increases, more steps are needed to integrate the trajectories 
accurately. Extensive tests of the code and comparisons with 
other numerical iV-body codes can be found in Kravtsov 
(1999) and Knebe et al. (2000). 

The code was modified to handle DE of different types, 
according to the prescription of Mainini et al. (2003b). Mod- 
ifications include effects due to the change in the rate of the 
expansion of the Universe and on initial conditions, keep- 
ing into account also spatial fluctuations of the scalar field 
before they enter the horizon. 

In this paper we use 4 new simulations. The mod- 
els are normalized assuming ug = 0.9. They are run in 
a box of 100/i _1 Mpc. We use 256 3 particles with mass 
m p — 4.971 • 10 9 hT 1 Mq. The nominal force resolution is 
3/i _1 kpc. All models are spatially flat, while f2 mo = 0.3 and 
h — 0.7. The two ACDM models start from different ran- 
dom numbers and are indicated as ACDM1 and ACDM2. 
The two DE models, named RP3 and SUGRA3, are started 
from the same random numbers of ACDM1. 



4 GALAXIES IN HALOS 

Halos made by more than 30 particles were found in sim- 
ulations by the same spherical overdensity (SO) algorithm 
used in Klypin et al (2003). The algorithm locates all non- 
overlapping largest spheres where the density contrast at- 
tains a given value A„. Density contrasts are assigned the 
virialization values, which depend on the redshift z and on 
parameters of the DE. For instance, at z = 0, A„ = 101.0 
for ACDM, 119.4 for SUGRA3 and 140.1 for RP3; A„ values 
for higher z can be found in Mainini, Maccio & Bonometto 
(2003; see also Mainini et al 2003b). 

Figure 3 shows the halo mass function in the LCDM1, 
the SUGRA3, and the RP3 models at z = 0. Differences be- 
tween the models are only due to different w(z), while their 




log (M/M h-») 

Figure 3. Halo mass function at z = 0. Results for 
ACDM (dashed line), RP (solid line), and SUGRA (dotted line) 
overlap. 

as is identical. Accordingly, at z = their mass functions 
almost overlap and are well fitted by a the Sheth-Tormen 
(ST, Sheth & Tormen, 1999) approximation. 

However, when we consider galaxies, we need to use 
another approach because individual halos may host many 
galaxies of different mass and luminosity. In order to assign 
galaxies to halos we use a HOD. This is a relatively novel 
approach to locate galaxies in each halo. It can be used in 
different ways. Full-scale semi-analytical methods can pre- 
dict such quantities as the luminosity, colors, star-formation 
rates. Unfortunately, many important mechanisms are still 
poorly understood making the results less reliable. It seems 
therefore advisable to minimize the physical input, keeping 
just to gravitational dynamics. 

In this paper we use a prescription consistent with the 
results of Kravtsov et al (2004). We utilize an analytical 
expression recently proposed by Vale & Ostriker (2004), but 
parameters of the approximation are different form Vale & 
Ostriker (2004) and tuned to produce a good fit to results 
of Kravtsov et al (2004). The approximation is based on 
the assumption that the probability P S (N S \M) for a halo of 
mass M to host N a subhalos is approximately universal. 
We take the Schechter approximation 

for the number of subhaloes with masses in the range m to 
rn + dm, for a parent halo of mass M. A must be such that 
the total mass in subhaloes, J"°° dmmN(m\M), is a fraction 
7M of the parent halo mass. Therefore A — j//3T(2 — a), 
so that the number of sub-halos of mass m is 

n sh {m)= I dM N{m\M)n h {M) , (7) 
Jo 

(n/i(M) is the halo mass function) independently from the 
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Figure 4. Comparison of different halo occupation distributions. 
The solid line is the HOD given by Kravtsov et al (2004). Other 
curves are obtained from eq.(8) for different parameters 7 as in- 
dicated in the plot. 



parent halo mass. The expression (6) yields the following 
number of subhaloes with mass > m in a halo of mass M: 



N sh (> m,M) = 



T f 

0T(2-«) j m 



dx x a exp(— x) . (8) 



/f3M 



Sub-haloes will be then identified with galaxies. Figure 
4 shows that the expression approximates the results of 
Kravtsov et al (2004) once we fit parameters a, f3, 7 and add 
to the expression (8) a unity, i.e. the halo as a sub-halo of 
itself. For large haloes, to be interpreted as galaxy clusters, 
this sub-halo could be the central cD; but adding an extra 
object, in such large galaxy sets, is just a marginal reset. For 
small haloes, instead, it is important not to forget that they 
represent a galaxy, as soon as they exceed the galaxy mass 
threshold. 

In a different context, Vale & Ostriker (2004) use the 
value 7 = 0.18, /3 = 0.39. Owing to the use we make of 
eq. (8), 7 = 0.7 appears more adequate. 

If the (differential) halo mass function rih(M) is known, 
the sub-halo mass function is 



N sh (> m) = 
/ dMn h (M) 

J m/7 



7 



(3T{2 - a) 

dxx~ a exp(— x) 

/0M 



(9) 



If we perform nonlinear predictions, n h (M) is obtained from 
the expression (5) or from the corresponding expressions in 
the ST case. When we deal with simulations, halo masses 
have discrete values m v = vm v , appearing times, up to 
a top mass VMfn v . Then 



N sh {> m) 



7 



A 



1 - 



1 1 1 1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 

- 


M* = 4xl0 IZ Mgh- 1 


\ • 
\ • 
\ • 
\ • 

\ • 


a= -0.8 


•s. - 


O = 1200 

, , , 1 , , , , 1 , , 





11 11.5 12 12.5 

Log 10 (M /M e h-') 



13 



Figure 5. Galaxy (cumulative) mass function at z = from 
simulations (dots) vs. an integral Schechter function, with <I> , a 
and M„ shown in the frame. 



dxx a exp(— x) 



(10) 



vm/f3M 



In Figure 5 we plot the galaxy mass function obtained 
with eq. (10), using the mass function of halos in the simu- 
lations, and identifying sub-haloes with galaxies. At z = 
model differences are unappreciable and the plotted function 
holds for all models. We also plot a Schechter function with 
the parameters shown in the frame, selected to minimize the 
ratios between differential values at all points. As expected, 
the two curves are close. In fact, there must be some rela- 
tion between masses and luminosities, but the M g /L g ratio 
should not be a constant. 



5 GALAXY ANGULAR DENSITY IN MODELS 
WITH DIFFERENT DE 

Let us now consider the galaxy mass function at higher z, 
for the different models. According to eq. (1), the number 
of galaxies with mass > m, in a solid angle A9 (A9 <JC 1), 
between redshifts z and 2 + Az, is 



iV g (> m,z;Az,A6) ^ r 2 (z) 
AzA6 2 ~ C H(z) 



n g (> m, z) 



(11) 



if n g (> m,z) is the comoving number density of galaxies 
with mass > m at a redshift z. The galaxy density can 
be obtained from the subhalo mass functions (9) and (10). 
Accordingly, the average angular distance 9 gg is given by 



e gg {>m,z)^A~z~ — 



H{z) 



n g (> m,z) 



1/2 



(12) 



/8T(2 - a) 



Thus, the l.h.s. expression is independent of the particular 
volume considered. 
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Such gg (> m,z) therefore depends on geometry, halo 
mass function and HOD. We however expect that the red- 
shift dependence mostly arises from geometry, while evolu- 
tion plays a significant role at higher z. In fact, the main 
difference between this and the cluster case is that evolu- 
tion is mild and discrepancies between models, in comoving 
volumes, up to z ~ 2, are even weaker. 

Let Oa{z), 6su(z), 0rp(z) be the mean angular dis- 
tances between galaxies at redshift z for the ACDM, SUGRA 
and RP models, respectively. Besides these functions, let 
us also consider the angular distance 9 geo (z) obtained from 
cq. (12), keeping the value of n g (> m, z = 0) at any redshift 
and the ACDM geometry. Therefore, 9 gso (z), although the 
symbol has no reference to ACDM, describes the behavior 
of the angular separation in a ACDM model, in the limit of 
no halo evolution. 

In Figure 6, we compare results obtained from simu- 
lations with the ST predictions in different models. Rather 
than presenting 9 mo d{z), we plot the fractional difference 
(#mod — 0geo)/0g eo . In principle, error bars can be evaluated 
in two ways: (i) by comparing ACDM1 with ACDM2 (cos- 
mic variance) and (ii) by comparing the differences between 
models. The latter can be done only at present because the 
models have the same power spectrum only at z = 0. The 
differences between models exist because the fluctuations 
grow differently in the past. At larger z, this evolutionary 
variance should be smaller, but is not easy to evaluate. We 
use differences between models at z — as a rough estimate 
of error bars at all redshifts. Judging by the differences be- 
tween ACDM1 and ACDM2, the cosmic variance seems to 
be smaller by a factor of three than the evolutionary differ- 
ences. We find similar behaviour for different galaxy masses. 
In all cases, differences between models can be clearly seen. 

The largest differences between models are attained at 
z ~ 1. Let us remind that the plot shows the fractional dif- 
ferences between DE models and the 2-dependence due to 
the mere ACDM geometry. At z ~ 1, this difference is just 
~ 5% for ACDM while, for SUGRA, it is ~ 20%, because 
of the different geometry and a still slower cosmological evo- 
lution. Even larger is the difference with RP. This compares 
with an evolutionary variance hardly exceeding ~ 2 %, if the 
effective comoving volume inspected is ~ 10 6 /i _3 Mpc 3 . For 
Az ~ 0.1, this corresponds to 89 ~ 30-40°. 

The discriminating power of this theoretical prediction 
is to be still compared with two possible sources of error: (i) 
Peculiar velocities, setting individual galaxies into an ap- 
parent redshift band different from the one to which they 
belong, (ii) Luminosity evolution. 

Overcoming the latter point is critical to the use of 
galaxies as indicators of DE nature and we shall devote the 
whole next section to the impact of luminosity evolution. 
The conclusion of this discussion is that galaxies are indeed 
a possible indicator of DE nature, but more work is needed 
before they can be efficiently used to this aim. 

In this section we shall report the result of a test per- 
formed to evaluate the impact of redshift displacements due 
to peculiar velocities. 

We divided the volume L 3 of the box into cells of side 
10 _1 L, with volume 10~ 3 L 3 . 50 such cells were selected at 
random at each redshift and replaced by the cells at the 
closest redshift considered with their galaxy contents. The 
redshift displacement between the original cell and its re- 
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Figure 6. The redshift dependence of the mean fractional angular 
separation in different models. Points with error bars show results 
of simulations. Curves indicate ST predictions. The plot shows a 
strong dependence of the galaxy redshift distribution on the DE 
nature. 



placements is ±0.5. For a change of 5 % in galaxy contents 
the average shift of 9 gg is ~ 1.2 %. 

This shift lays well below the error bars shown in Figure 
6. However, the average errors due to evolutionary and cos- 
mic variance plus peculiar velocities are ~ 2.4 % and the 
whole error never exceeds 2.6-2.7 % . We therefore argue 
that these sources of error do not affect the robustness of 
results. 



6 DEPENDENCE ON MASS-LUMINOSITY 
RELATION 

Let us now consider galaxy evolution, which is expected to 
cause 2-dependence of the average M g /L g ratio, but can 
also yield a 2-dependence of the M g /L g distribution about 
such average, in a way which may depend on the mass range 
considered. From an observational point of view, when we 
consider galaxies of various luminosities L g , we must then 
take into account that their expected mass M g could be 
distributed with different laws at different z and L g . 

Let us therefore consider the galaxy distribution on the 
M g , Lg plane at a given z, yielding the galaxy number 

dN = D{Mg,Lg,[z])dMgdLg (13) 

in the infinitesimal area dM g dL g about the point M g , L g . 
We put 2 in brackets to outline that, in respect to it, D 
is not a distribution but a function . Obviously we expect a 
strong correlation between M g and L g , at any z, so that it 
makes sense to consider an average M g /L g ratio. 

Once the distribution D is assigned, the distributions 
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on M g ( at fixed L g ) and on L g (at fixed M g ) read 
4>{M„,[z])= [ dL g D(M g ,L g ,[z]) , 



^L g ,[z]) 



J dM g 



D(M g ,L g ,[z]). 



The number dN is also the product of <p(M g , [z]) times the 
distribution on luminosities at fixed mass M g : 

dN = 4>{M g , [z]) Q(L g ; [M g , z]) dM g dL g . (14) 

Equating the r.h.s.'s of eqs. (13) and (14) yields 

D(M g ,L g ,[z]) 



Q(L g ;[M g ,z]) 



(15) 



/ dlD(M g ,l,[z}) 
and, similarly, the distribution on masses at given L g , reads 
D(M g ,L g ,[z]) 



P(M g ;[L g ,z}) 



J dm, D(m, L g , [z]) 



(16) 



We can now use P to work out the average M g /L g at fixed 
L g , and the distribution on M g /L g about such average. 
Clearly 



Ma 



1 



dmmP(m; [L g , z]) = 



(17) 



(18) 



1 J dram D(m, L g ,[z\) 
L g J dm D(m, Lg,[z]) 
while the distribution 

V(M B ;[L„z]) = ¥L D(M 9 ,L B ,[z]) 

L g J dmD(m, L g , \z\) 

tells us how M g /L g is distributed around (M g /L g ) L z . 

The impact of the evolution of stellar populations (or 
other mechanisms) on the M g /L g ratio can be fully ex- 
pressed through the distribution D(M g ,L g ; [z]) in eq. (13). 
From it we can work out an average mass/luminosity ra- 
tio {Mg/Lg} and the distribution on masses V; they both 
depend on L g and z. 

Let us now try to inspect how such variable T> distribu- 
tion affects our results, taking into account that we mostly 
ignore how such variations occur. Accordingly, we shall pro- 
ceed as follows: we define a "wild" distribution, that we ex- 
pect to spread the M g /L g ratio, at fixed L g , farther from 
average than any physical V, at any z, will do. The effects 
caused by such wild distribution should then be an overes- 
timate of the effects of the actual distributions. Should they 
cause just a minor perturbation in estimates, all we have 
to worry about is the redshift dependence of the average 

(Mg/Lg) Lg . 

More in detail, we shall allow that a galaxy of given 
mass M g has a luminosity in an interval L\,L^ with L2 ~ 
20 L\ . We test this prescription without direct reference to 
luminosities: In the sample of galaxies obtained through the 
HOD, at each z, each galaxy mass (m) is replaced by a mass 
m' — m + AmR, R being a random number with normal 
distribution and unit variance. We take Am = 0.8 m, but 
replace all m! < 0.1 m with 0.1 m, as well as all m' > 1.9 m 
with 1.9 m. The shift is therefore symmetric on m (not on 
logm). The operation causes a slight increase of the mass 
function above M ~ 2.8 • 1Q 11 h -1 M@ (by a few percents), 



as there are however more lighter galaxies coming upwards 
than heavier galaxies going downwards (below M , the low- 
mass cut-off of the mass function, set by the mass resolution 
of our simulations - see section 3 - begins to cause shortage 
of transfers upwards). The operation is then completed by 
reducing all masses by a (small) constant factor, so that, 
summing up all masses of objects with mass M g > M, we 
have the same total mass as before the operation. This lowers 
the limit below which the mass function preserves its initial 
shape, but we never use galaxy samples including masses 
below 3- W 11 h- 1 M e . 

We re-estimated 8 mod \/ Az using the new masses, for 
the same mass limits as before, and compare the changes 
obtained in this way with the Poisson uncertainty, due to 
the finite number of galaxies in each sample. 

We find that the error obtained from the above proce- 
dure ranges between 20 and 40 % of the Poisson error. 

This output tells us that the evolution of the physi- 
cal distribution can be expected to redistribute results well 
inside Poisson uncertainty. If we expect this redistribution 
to be random, the top value of the whole expected error is 
then 3 %; if we attribute a systematic character to it and re- 
frain from performing a quadratic sum with the other error 
sources considered, the overall possible error is still within 
3.8 %. It should be outlined that here we pushed all error 
sources to their maximum; thus, we believe that the above 
estimates are safely conservative. 

Let us now discuss how the evolution of the average 
M g /L g ratio can affect the use of galaxies to detect DE 
nature. In principle, one could use the z dependence of 
D(M g , L g ), to obtain the z-dependence of the (M g /L g )L g 
ratio. More realistically, suitable data sets can directly pro- 
vide the z-dependence of (M g / L g ) l„, with some residual 
uncertainty. The basic issue is then: How well the z evolu- 
tion of {M g /L g ) L g is to be known, in order that we can test 
DE nature ? 

Figure 6 is however devised so to provide a direct reply 
to this question: If we double the values of 56 g g/6 gg pro- 
vided there, we have a fair estimate the difference between 
evolution rates of TZ = (M g /L g ) Lg (z)/(M g /L g ) Lg (z = 0) 
needed just to cover the differences between geometry and 
dynamics for the ACDM model or to compensate the differ- 
ences between the ACDM geometry and whole evolution for 
dynamical DE models, (in the limit M g <C M», where M* 
is the mass scale appearing in a Schechter-like expression). 
Figure 6 can just be interpreted in both senses, changing the 
name of the ordinate. 

For instance, at z = 0.5, an uncertainty ~ 20 % (35 %) 
on the evolution of M g /L g is needed to hide the difference 
between ACDM and SUGRA (RP). 



Let us show this point. Owing to eq. (12), 6 gg oc n g 



1/2 



so that a shift 80 gg in the observed angular distance arises 
from a shift 



8n„ 



JO 99 



in the galaxy number density 



M> Mg) =n g ( > L g ^- 



(19) 



(20) 



The latter shift, as shown in eq. (20), can arise from a shift 
on Mg/Lg, being 
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8n g = n g (M g )L g 5 



Ma 



(21) 



here n g (M g ) is the differential mass function, obtained by 
differentiating the integral mass function n g (> M g ). There- 
fore, 



Sn g (>M g ) ^ S(M g /L g ) 
M g n g {M g ) (M g /L g ) 

and 

S(M g /L g ) 8n g (> M g ) n g (> M g ) 



(22) 



(23) 



{M g /L g ) n g (> M g ) M g n g {M g ) ' 

If we approximate the integral mass function by a Schechter 
expression, it is \mn(m)/n(> m)\ = 1 + m/M t , so that 

8(M g /L g ) _ Sn g (>M g ) f M a r' 
(Mg/L g ) n g {>M g ) \ + Mj • 1 ' 

Using this equation together with eqs. (19) and (23), we 
have the relation 



S{M 9 /L g ) 

(Mg/Lg) 



of, , M g y 1 89 gg 
\ Mj 9 qg 



(25) 



telling us how to use Figure 6 to estimate the evolution of 
M g /L g needed to yield the same effects of a change in DE 
nature. This equation also tells us how to use Figure 6 for 
masses approaching M„ . 



7 THE VOID PROBABILITY FUNCTION 

Let us randomly throw spheres of radius R, in a space where 
objects of various masses M are set. The probability P (R) 
of finding no object with M > M tr in them, if the VPF for 
objects of mass > M tr . 

We expect and find no model dependence in the galaxy 
VPF's at z — 0. At z > 0, a critical issue is how M tr is set. 
One can simply plan to determine the galaxy masses M g 
from data (e.g., from L g values), so to select galaxies with 
M g > M tr - As widely outlined, this choice involves several 
complications. Another option is taking the most luminous 
galaxies up to an average angular distance 6 gg . 

Each threshold M tr , for any z and Az, yields a value 
of 9 gg . Fig. 6 shows how 9 gg depends on the model at fixed 
threshold. Viceversa, if we keep, for that z and Az, a fixed 
9 gg , the relative M tr varies with models. We can compare 
models either at fixed M tr or at fixed 9 gg . Dealing with 
observations, the latter option is easier, but mixes up the 
intrinsical VPF dependence on the model and other features, 
which also depend on the model. 

Besides of threshold setting, another issue bears a great 
operational relevance. In principle, VPF's can be evaluated 
in the comoving volumes where galaxies are set and com- 
pared there with VPF's from data. This is however unade- 
quate to evaluate how discriminatory is the VPF statistics. 
To do so, we rather follow the following steps: 

(i) From cartesian coordinates Xi (i = 1,2,3) in comoving 
volumes we work out the redshift and the celestial coordi- 
nates z, 9, <f> that an observer, set at z = 0, would measure. 
This is done by using the geometry of the model. 

(ii) Data also give z, 9, <f>, for each galaxy. But, to estimate 
the VPF, they must be translated into cartesian coordinates 
Xi. An observer can only perform such translation by using 



the geometry of a fiducial model, e.g., ACDM. The second 
step, to forge predictions, therefore amounts to re-transform 
z, 9, <f> into cartesian coordinates Xi, but using now the fidu- 
cial ACDM geometry; a;»'s coincide with the Xi's only for a 
ACDM cosmology. Let us call fiducial space, the environ- 
ment where galaxies are now set. 

(iii) We then estimate the VPF's, for all models, in the fidu- 
cial space; these VPF's should be compared with observa- 
tional data, but can also be compared one another to assess 
how discriminatory this statistics can be. 

Although comparing predictions for VPF's in comoving 
volumes, therefore, bears little discriminatory meaning, our 
outputs are more easily explained if we start from comoving 
space VPF's. Let us recall that, on galaxy scales, evolution- 
ary differences, up to z ~ 2, are modest. Accordingly, for 
an assigned M tr , we find just marginal discrepancies, as is 
shown in Fig. 7 for M tr = 6- 10 11 M^k -1 . Differences among 
VPF's arise, of course, if we do not fix Mtr, but 9 gg , as a 
reflex of 9 gg differences. These VPF's are shown in Fig. 8. 

VPF's in respect of comoving coordinates bear a strict 
analogy with mass functions in comoving volumes. The fact 
that geometry erases almost any signal on the cluster mass 
function is analogous to what happens for the VPF, when 
we pass from the comoving to the fiducial space. This fact is 
far from trivial. Let us compare these VPF's with the VPF's 
in a Poisson sample with the same 9 gg (Fig. 10) and remind 
that (White 1979) 



Po(R) = exp 



(-N R r An) 



e n> (R) 



(26) 



Here Nr is the average number of points in a sphere of radius 
R; ^ n \R) are the n-point functions averaged within the 
same sphere. For the Poisson sample P (R) = exp(— Nr), 
as all ^ n \R) vanish. The difference between Poisson VPF 
and model VPF is to be fully ascribed to , as Nr is set 
equal. This difference is huge, in respect to the differences 
between models, which should arise because of shifts. 
Their paucity indicates that density renormalization almost 
erases the shifts in correlation functions of all orders. 

The cancellation between geometrical and 9 gg effects, 
shown in Fig. 9, indicates that the passage from comoving 
to fiducial coordinates bears a weight comparable with the 
differences shown in Fig. 8. It therefore comes as no surprise 
that VPF's, for fixed Mtr, almost absent in the comoving 
space, are significant in the fiducial space. They are shown 
in Fig. 11 and, as is expected, the curves of different models 
appear in the opposite order in respect to Fig. 8. 

If the redshift dependence of the M g /L g ratio is un- 
der control, Fig. 11 shows a discriminatory prediction which 
can be compared with data. Once again, the problem con- 
cerns both the evolution of the mean M g /L g ratio and sin- 
gle galaxy deviations from average. Supposing that the av- 
erage M g /L g evolution is under control, we can estimate 
the impact of individual deviations by replacing the sharp 
threshold on M g with a soft threshold, substituting each 
galaxy mass M g with M g + AM g R, as in the previous sec- 
tion. This test was performed for samples as wide as those 
in our 100/i _1 Mpc box and the effects of such replacement 
are modest, amounting to ~ 10% of the difference found 
between ACDM and SUGRA. 

Accordingly, the critical issue concerns the mean M g /L g 
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Figure 7. VPF's in comoving volumes for M tr = 6- W 11 M & h~ l . 
The four panels refer to redshift 0.5, 1, 1.5, 2, as indicated in the 
frames. The galaxy numbers in the simulation box are reported, 
for each model and redshift. Solid, dashed and dotted lines as in 
Fig. 6. 
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Figure 8. VPF's in comoving volumes for fixed angular density. 
The four panels refer to redshift 0.5, 1, 1.5, 2, as indicated in the 
frames. Numbers in the frames as in Fig. 7. Solid, dashed and 
dotted lines as in Fig. 6. 



ratio. By comparing VPF's for different thresholds, we can 
however see that, in order that the shift of M tr induces a 
VPF shift similar to differences between models, Mtr is to 
be displaced by a factor 1.7-1.8. Uncertainties ~ 10% on 
the mean M g /L g ratio would therefore leave intact the dis- 
criminatory power of the VPF statistics. 

Before concluding this section, let us finally comment 
on sample variance. All ACDM VPF's were estimated on 
the ACDM1 simulation. Differences between ACDM1 and 
ACDM2 are however small and could not be appreciated 
in the above plots. Accordingly, sample variance is not a 
relevant limit to the use of VPF's. 



8 CONCLUSIONS 

Dark Energy modifies the rate of cosmic expansion in the 
epoch when a substantial fraction of fluctuations on clus- 
ter scales reach their turnaround. Therefore, it seems quite 
natural to trace the redshift dependence of w(z) = pde/pde 
using the cluster mass function at different redshifts. 

Unfortunately, the situation is more complicated. The 
evolution of w(z) affects observations in two ways. First, it 
causes objects to form and evolve with different rate. Second, 
it results in a different mapping of comoving coordinates of 
galaxies to observed angular positions and redshifts. Then, 
on a cluster mass scale (~ 10 14 /i _1 M Q ), the evolutionary 
and geometrical effects tend to cancel, which makes clusters 
somewhat problematic for testing the equation of state. 

In this paper we discuss the use of scales where the 
evolutionary effects are minimized, so that the geometrical 
effects leave a clearer imprint. In a sense, this is not a new 
procedure: an analogous idea is utilized when the DE equa- 
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Figure 9. VPF's in fiducial volumes for fixed angular density. 
Numbers in the frames as for Fig. 7. Solid, dashed and dotted 
lines as in Fig. 6. 

tion of state is tested by using a standard candle. Obviously, 
galaxies are not standard candles themselves, but they can 
provide a "standard meter" through their (almost) model- 
independent abundance and evolution. 

Previous analysis, which focused on clusters, tried to 
overcome the above difficulties by making recourse to vari- 
ous features. For example, the evolutionary dependence on w 
is preserved if one considers masses weM above 10 14 /i _1 M Q . 
Unfortunately, clusters with masses ~ 10 15 ft _1 M Q or lar ger 
are rare today and surely are even more rare in the past. 
Some analysis (see, e.g., Haiman et al. 2000) stressed a pos- 
sible role of very massive clusters at large z. Yet, the number 
of such clusters cannot be large and, thus, comparing such 
predictions with observations is a hard challenge. There is 
also another problem with using cluster masses. Usually the 
mass function is estimated with a PS-like approximations, 
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Figure 10. Differences between VPF's of various models and 
the VPF for a Poisson sample, in fiducial volumes for fixed an- 
gular density. They arise from the sum of n-point correlation 
functions in eq. (26). Tiny residual differences between models, 
almost unappreciable in the previous plot, arise from differences 
between their n-point functions, clearly almost erasen by geomet- 
rical rcnormalization. 
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Figure 11. VPF's in fiducial volumes for fixed mass limit (see 
text). Numbers in the frames as for Fig. 7. Solid, dashed and 
dotted lines as in Fig. 6. 



which are well tested with simulations. The "virial" radius 
R v is defined so that inside R v the density contrast is A„. 
The value of A„ depends on the redshift and on DE model. 
On the contrary, data are typically analized with a standard 
density contrast A c ~ 180 (or 200). Increasing A c reduces 
the amplitude of the mass function. If mass functions de- 
fined with variable A„ (almost) overlap one another, mass 
functions defined with constant A c can be different. In order 
to account for these differences in the definitions, one needs 
to assume some shape for the density profile in the outskirts 
of clusters. This is typically done by using a NFW profile 
with concentration c s ~ 5. As we deal with rather peripheral 
(virial) cluster regions, we can neglect the spread of actual 
values of concentration. However, when different w(z) are 



considered, the mean c a changes substantially, up to 80% 
(Klypin et al 2003; Kuhlen et al. 2005). The differences be- 
tween M200 and M v are not large - 10-15 %. Approximately 
the same percent of the difference depends on w. This may 
be still important. Neglecting these corrections may lead to 
substantial systematic errors. 

As an alternative to using galaxy clusters, so avoiding 
these and other problems, here we suggest to exploit the de- 
pendence on DE nature of the redshift distribution of galax- 
ies and argue that the difficulties of this approach can be 
overcome. 

A first problem is that one needs to know how to treat 
subhalos of more massive halos because a large fraction of 
galaxies is hosted by subhalos. To populate massive halos 
with galaxies we use recent results on the halo occupation 
distribution (HOD). Accordingly, we believe that this diffi- 
culty can be readily overcome. 

Dealing with galaxies also requires knowledge of their 
masses M a . The M g — L g relation is more complex than the 
relation between M v , X-ray flux and T in clusters. 

Galactic evolution studies and/or techniques aiming to 
compare dynamical or lensing masses with luminosities can 
be used to this aim (Bressan et al 1994, Portinari et al 2004, 
Prada et al 2003). It is also known that there is no one-to- 
one correspondence between L g and M g , as the luminosities 
of two galaxies of the same mass can be quite different, up to 
one order of magnitude. We then considered separately two 
different issues: (i) How well we must know the evolution 
of (Mg/Lg)L g - (u) Which can be the impact of fluctuations 
about such average value, taking into account that the dis- 
tribution about average can depend on z and mass. 

How precisely (M g /L g ) L g is to be known, in order that 
different cosmologies can be safely discriminated, is shown 
by Fig. 6, according to Sec. 6: If we double the values of 
59 gg /8 gg provided there, we have a fair estimate the evolu- 
tion rate TZ = (M g / L g ) Lg (z) / (M g / L g ) Lg (z = 0) just cov- 
ering the differences between models. Typically, if the es- 
timated {M g I 'L g ) L g evolution is reliable at a ~ 10% level, 
different models could be discriminated. 

To reach this goal, fresh observational material is 
needed, but no conceptual difficulty is apparently involved 
in its acquisition. 

The spread of the M g /L g ratio, around its average 
value, also involves a delicate issue, as it could cause sys- 
tematics. Here we reported the effect of a random spread of 
L g in a luminosity interval L\ , L2 with L2 — 20 L\ . If the 
physical distribution of luminosities, for any mass and at any 
redshift, is within these limits, then possible systematics are 
well within errors due to other effects. 

Further tests on L g spread were also performed, which 
are not reported in detail in this paper. They apparently 
indicate that really wide and ad-hoc distributions are nec- 
essary, in order that possible systematics exceed Poisson un- 
certainty. 

Although bearing these reserves in mind, we conclude 
that estimates of the redshift dependence of the average 
Mg/Lg, reliable within ~ 10%, can enable us to obtain a 
fair information on DE nature. 

In this paper we also discuss various tests based on the 
void probability function. We find that the VPF is almost 
model independent when estimated for samples with con- 
stant angular number density of galaxies. This result ap- 
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parently suggests that n-point functions are almost model 
independent, at any z, once distances are suitably rescaled. 
If the mean M g /L g evolution is under control, at the 10% 
level, we can also assess that VPF, for samples defined with a 
given Mtr, puts in evidence geometrical differences between 
models and provides a discriminatory statistics. 

Combining the simulated halo distribution with the 
HOD provides an effective tool for testing the equation of 
state of DE. More work is needed to define the L g -M g rela- 
tion and, possibly, to reduce systematic effects. It is justified 
to expect that the 2-dependence of galaxy distribution, in 
deep galaxy samples, will allow us to constrain the DE na- 
ture even more reliably than the density of galaxy clusters 
in future complilations, sampling them up to large z's. 
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APPENDIX A: DYNAMICAL DARK ENERGY 
MODELS 

Dynamical DE is to be ascribed to a scalar field, <j>, self- 
interacting through an effective potential V(<f>), whose dy- 
namics is set by the Lagrangian density: 

Cde = -\^9 (d^d^ + V(4>)) . (Al) 

Here g is the determinant of the metric tensor g M „ = 
(^(r^jdx^dxv (r is the conformal time). In this work we 
need to consider just a spatially homogeneous (f> (di4> <C <f>; 
i = 1,2, 3; dots denote differentiation with respect to r); the 
equation of motion is then: 

+ 2-0 + a 2 ^=O. (A2) 
a d(p 

Energy density and pressure, obtained from the energy- 
momentum tensor T M „, are: 

P=-T$ = ^+V{4>), p=It? = H-V(0, (A3) 
so that the state parameter 

«, s g = £^zw (A4) 

P 4> 2 /2a + V(<t>) 

changes with time and is negative as soon as the potential 
term V(<j>) takes large enough values. 

The evolution of dynamical DE depends on details of 
the effective potential V(<j>). Here we use the model proposed 
by Ratra & Peebles (1988), that yield a rather slow evolution 
of w, and the model based on supergravity (Brax & Martin 
1999, 2001; Brax, Martin & Riazuelo 2000). The latter gives 
a much faster evolving w. The RP and SUGRA potentials 

A 4+a 

V W = — RP > ( A5 ) 
A 4+c * 

V{4>) = —— exp(47rG<£ 2 ) SUGRA (A6) 



cover a large spectrum of evolving w. These potentials allow 
tracker solutions, yielding the same low-z behavior that is 
almost independent of initial conditions. In eqs. (A5) and 
(A6), A is an energy scale in the range 10 2 -10 10 GeV, rel- 
evant for the physics of fundamental interactions. The po- 
tentials depend also on the exponent a. Fixing A and a, the 
DE density parameter Qde,o is determined. Here we rather 
use A and O,de,o as independent parameters. In particular, 
numerical results are given for A = 10 3 GcV. 

The RP model with such A value is in slight disagree- 
ment with low-Z multipoles of the CMB anisotropy spectrum 
data. Agreement may be recovered with smaller A's, which 
however loose significance in particle physics. The SUGRA 
model considered here, on the contrary, is in fair agreement 
with all available data. 
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